function Q_new = findQ2_v2(Q, i, j, P, n)
    Q_new = Q;
    Q1 = Q(i,:) * (1-P(j));
    Q2 = [0, Q(i,1:n)] * P(j);
    Q_new(i, :) = Q1 + Q2;
return;